Nodal superconductivity in Ba(Fei_ 2) Ru x )2As 2 induced by isovalent Ru substitution 
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We present the ultra-low-temperature heat transport study of an iron-based superconductor 
Ba(Fe .64Ruo.36)2As 2 (T c = 20.2 K), in which the superconductivity is induced by isovalent Ru 
substitution. In zero field we find a large residual linear term Ko/T, more than 40% of the normal- 
state value. At low field, the ko/T shows an H 1 ^ 2 dependence. These provide strong evidences for 
nodes in the superconducting gap of Ba(Fe .64Ruo.36)2As2, which mimics that in another isovalently 
substituted superconductor BaFe2(Asi- iE P a ;)2. Our results show that the isovalent Ru substitution 
can also induce nodal superconductivity in BaFe2As2, as P does, and they may have the same origin. 
We further compare them with other two nodal superconductors LaFePO and LiFeP. 

PACS numbers: 74.70.Xa, 74.25. fc, 74.20. Rp 
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Since the discovery of high-T c superconductivity in 
iron-based compounds [l|, 0] , the electronic pairing mech- 
anism has been a central issue 0] . One key to understand 
it is to clarify the symmetry and structure of the super- 
conducting gap 0]. However, even for the most studied 
(Ba,Sr,Ca,Eu)Fe2As2 (122) system, the situation is still 
fairly complex [J]. 

Near optimal doping, for both hole- and electron- 
doped 122 compounds, the angle-resolved photon emis- 
sion spectroscopy (ARPES) experiments clearly demon- 
strated multiple nodeless superconducting gaps @, 0, 
which was further supported by bulk measurements such 
as thermal conductivity ■ On the overdoped side, 
nodal superconductivity was found in the extremely hole- 
doped KFe2As2 HJ U|i while strongly anisotropic gap 
or isotropic gaps with significantly different magni- 
tude [Hi Gil were suggested in the heavily electron-doped 
Ba(Fei_ :c Co a ;)2As2. On the underdoped side, recent heat 
transport measurements claimed possible nodes in the su- 
perconducting gap of hole-doped Bai_ a; K 2 ;Fe2As2 with 
x < 0.16 [l4j| . in contrast to the nodeless gaps found in 
electron-doped Ba(Fei_a; Co^ As2 [§]. 

Intriguingly, nodal superconductivity was also found 
in BaFe 2 (Asi_ K P 2; )2 (T c = 30 K) [lH [lj], in which 
the superconductivity is induced by the isovalent P sub- 
stitution for As. The laser-ARPES experiments on 
BaFe2(Aso.65Po.35)2 showed isotropic gaps in the three 
hole pockets around the Brillouin zone (BZ) center, 
therefore the gap nodes must locate on the electron pock- 
ets around the BZ corners 17J. Moreover, previously 



LaFePO (T c ~ 6 K) displays clear nodal behavior [18- 
2fij . and recently there is penetration depth evidence for 
nodes in the superconducting gap of LiFeP (T c ~ 4.5 K) 
2l| . The nodal superconductivity in these P-substituted 
compounds are very striking, which raises the puzzling 
question why the P substitution is so special in iron- 
based superconductors. The theoretical explanations of 
this puzzle are far from consensus 22-25)]. 



In the Fe2As2 slabs of iron arsenides shown in Fig. 1, 
instead of substituting As with P, there is an alternative 
way for isovalent substitution, to substitute Fe with Ru. 
Indeed, superconductivity with T c up to 20 K was found 
in Ba(Fei_ 2; Ru :E )2As2 and Sr(Fei_ a; Ru :!; )2As2 [13]. 
The phase diagram of Ba(Fei _ a Rug )oAs2 is very simi- 
lar to that of BaFe 2 (Asi_ :i; P a; )2 [H, The ARPES 
measurements on Ba(Feo.65Ruo.3s)2As2 showed that Ru 
induces neither hole nor electron doping, but the hole and 
electron pockets are about twice larger than in BaFe2As2 
[30l ] . To investigate the superconducting gap structure in 
these Ru-substituted superconductors may help to solve 
above puzzle of P substitution. 

In this Letter, we report the demonstration of 
nodal superconductivity in optimally substituted 
Ba(Feo.64Ruo.36)2As2 by thermal conductivity mea- 
surements down to 50 mK. Our finding shows that 
the nodal superconducting states in P-substituted iron 
arsenides are not that special, and suggests a common 
origin of the nodal superconductivity induced by isova- 
lent substitutions, at least in Ba(Fei_ a; Ru a ;)2As2 and 
BaFe2(As!_ x P x )2. 

Single crystals of optimally substituted 
Ba(Feo.64Ruo.36)2As2 were grown using a self-flux 
method [31(. Plate-shaped crystals with shiny surfaces 
were extracted mechanically. The Ru substituting 
level was determined by energy dispersive X-ray spec- 
troscopy. The dc magnetic susceptibility was measured 
at H = 10 Oe, with zero-field cooled, using a SQUID 
(MPMS, Quantum Design). The sample was cleaved 
to a rectangular shape of dimensions 1.50 x 0.68 
mm 2 in the afr-plane, with 70 /xm thickness along the 
c-axis. Contacts were made directly on the sample 
surfaces with silver paint, which were used for both 
resistivity and thermal conductivity measurements. The 
contacts are metallic with typical resistance 200 mil 
at 1.5 K. In-plane thermal conductivity was measured 
in a dilution refrigerator, using a standard four-wire 
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FIG. 1: (Color online). Crystal structure of BaFe2As2. There 
are two ways for isovalent substitution in the Fe2As2 slabs, 
substituting As with P, or Fe with Ru. Both substitutions 
can induce superconductivity, and result in similar phase di- 
agrams. 



steady-state method with two Ru02 chip thermometers, 
calibrated in situ against a reference Ru02 thermometer. 
Magnetic fields were applied along the c-axis. To ensure 
a homogeneous field distribution in the sample, all fields 
were applied at temperature above T c . 

Fig. 2a shows the in-plane resistivity p(T) of our 
Ba(Fe .64Ruo.36)2As2 single crystal. The magnitude and 
shape of p(T) are consistent with previous report [28j |. 
The normalized magnetization was plotted in the in- 
set, which displays a nice superconducting transition 
at about 20 K. According to the phase diagram of 
Ba(Fe 1 _ :c Ru :c ) 2 As2 [HHH, no static magnetic order ex- 
ists in our optimally substituted sample. The resistivity 
data between 30 and 90 K are fitted to p{T) = p + AT n , 
which gives a residual resistivity po — 43.7 ± 0.1 pficm 
and n — 1.31 ± 0.1. Such a non-Fermi-liquid temper- 
ature dependence of p(T) is similar to that observed in 
BaFe 2 (Asi_ a; P 2: )2 near optimal substitution, which may 
reflet the presence of antiferromagnetic spin fluctuations 
near a quantum critical point [32|. 

The low-temperature part of p(T) is plotted in Fig. 
2b. The zero-resistance point of the resistive tran- 
sition is at T c = 20.2 K, which is in good agree- 
ment with the diamagnetic superconducting transition 
shown in the inset of Fig. 2a. To estimate the up- 
per critical field H c2 , the resistivity in H — 3, 6, 9, 
and 14.5 T with H\\c were also measured, as seen in 
Fig. 2b. Fig. 2c shows the temperature dependence 
of H C 2, defined by p — and p — Q.lpN, respec- 
tively. For comparison, the H C 2 (T) of the electron-doped 
Ba(Feo.gCoo.i)2As 2 sing le crystal with T c « 22 K is re- 
produced from ref. |33j. One can see that although the 
T c of our Ba(Feo.64Ruo.36)2As 2 is only slightly lower than 
Ba(Feo.gCoo.i)2As 2 , its H C 2(T) is significantly lower. For 
Ba(Feo.gCoo.i)2As 2 , extrapolation of the H C 2(T) data 
suggests H C 2(0) between 40 and 50 T, much larger than 
that obtained from Werthamer-Helfand-Hohenberg for- 
mula H^ HH {0) = -0.69T c dH c2 /dT\ T = Tc Wk- For our 
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FIG. 2: (Color online). (a) In-plane resistivity of 
Ba(Feo.64Ruo.36)2As2 single crystal. The solid line is a fit 
of the data between 30 and 90 K to p(T) = po + AT n . Inset: 
normalized magnetization, (b) Low-temperature resistivity 
in H = 0, 3, 6, 9, and 14.5 T with H\\c. (c) Temperature 
dependence of the upper critical field H C 2- The squares and 
circles represent H C 2 of Ba(Feo.64Ruo.36)2As2 defined by p — 
and p — O.lpjv, respectively. The triangles represent H C 2 of 
Ba(Fe .9Coo.i)2As 2 with T c w 22 K [11]. 



Ba(Fe .64Ruo.36)2As 2 , H^ HH (0) 15.3 T is obtained, 
which is also apparently lower than the actual -ff C 2(0). 
We can only roughly estimate the bulk i/ c2 (0) « 23 T, 
defined by p — 0, by linearly extrapolating the data be- 
tween 6 and 14.5 T in Fig. 2c. Note that a slightly 
different H C 2(0) does not affect our discussion below. 

Fig. 3 shows the temperature dependence of the in- 
plane thermal conductivity for Ba(Feo.64Ruo.36)2As2 in 
H = 0, 1, 2, 4, 6, 9, and 12 T magnetic fields, plotted 
as k/T vs T. All the curves are roughly linear, as pre- 
viously observed in BaFei.gNi .i As2 0, KFe2As2 ic|. 
and overdoped Ba(Fei_a; Co^ As2 single crystals 0, 12j j. 
Therefore we fit all the curves to k/T = a + bT a ~ 1 with 
a fixed to 2. The two terms aT and bT a represent contri- 
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FIG. 3: (Color online). Low-temperature in-plane thermal 
conductivity of Ba(Feo.64Ruo.3e)2 As2 in magnetic fields ap- 
plied along the c-axis (H = 0, 1, 2, 4, 6, 9 and 12 T). The 
solid lines are n/T — a + bT fit to all the curves, respectively. 
The dash line is the normal-state Wiedemann-Franz law ex- 
pectation Lo/po, with Lo the Lorenz number 2.45 x 10~ 8 
WftK" 2 and p = 43.7 pQcm. 



butions from electrons and phonons, respectively. Here 
we only focus on the electronic term. 

In zero field, the fitting gives a residual linear term 



K /T 



0.266 ± 0.002 mW KT 2 cm" 



This value 



is more than 40% of the normal-state Wiedemann- 
Franz law expectation kno/T = L /p = 0.56 mW K -2 
cm -1 , with L the Lorenz number 2.45 x 10" 8 Wf2K~ 2 
and po — 43.7 /if2cm. For another isovalently substi- 
tuted BaFe2(Aso.67Po.33)2 single crystal, similar value of 
k /T ss 0.25 mW K~ 2 cm -1 was obtained, which is about 
30% of its normal-state kno/T [Hf. The significant kq/T 
of Ba(Fe .64R.uo.36)2As2 in zero field is a strong evidence 
for nodes in the superconducting gap [34 1 . 

The field dependence of kq/T may provide further 
support for the nodes 34]. In Fig. 4, the normal- 
ized (kq/T)/(kno/T) of Ba(Feo.64Ruo.36)2As 2 is plotted 
as a function of H/H C 2, with the normal-state knq/T 
= 0.56 mW K~ 2 cm" 1 and bulk H c2 = 23 T. Simi- 
lar data of the clean s-wave superconductor Nb (35} , an 
overdoped d-wave cuprate superconductor Tl-2201 [36| . 
and BaFe2(Aso.67Po.33)2 [HI are also plotted for com- 
parison. For a nodal superconductor in magnetic field, 
delocalized states exist out the vortex cores and dom- 
inate the heat transport in the vortex state, in con- 
trast to the s-wave superconductor. At low field, the 
Doppler shift due to supcrfluid flow around the vortices 
will yield an H 1 / 2 growth in quasiparticle density of 
states (the Volovik effect |37[), thus the H 1 ^ 2 field depen- 
dence of kq/T. From Fig. 4, the behavior of kq(H)/T in 
Ba(Feo.64Ruo.36)2As2 clearly mimics that in Tl-2201 and 
BaFe 2 (As .67Po.33)2. In the inset of Fig. 4, the k (H)/T 
of Ba(Feo.64Ruo.36)2As2 obeys the H 1 / 2 dependence at 
low field, which supports the existence of nodes in the 



FIG. 4: (Color online). Normalized residual linear term 
ko/T of Ba(Feo.64Ruo.36)2As2 as a function of H/H C 2- Sim- 
ilar data of the clean s-wave superconductor Nb J35| , an 
overdoped d-wave cuprate superconductor Tl-2201 [3g|, and 
BaFe2(Aso.67Po.33)2 [Rj are also shown for comparison. The 
behavior of k,q(H)/T in Ba(Feo.64Ruo.36)2As2 clearly mimics 
that in Tl-2201 and BaFe2(Aso.67Po.33)2- Inset: the same 
data of Ba(Feo.64Ruo.36)2As2 and Tl-2201 plotted against 
(H I H C 2) 1 '' 2 . The lines represent the H 1 ^ 2 dependence. 



superconducting gap. 

To our knowledge, previously there are five iron- 
based superconductors displaying nodal superconductiv- 
ity, KFe 2 As 2 [13, Hlj . underdoped Bai_ a; K :E Fe 2 As 2 (x < 
0.16) [3], BaFe 2 (Asi_ :c P a; )2 [H, (3, LaFePO ^mM, 
and LiFeP [2l| . Here we only consider the "in-plane 
nodes", not counting the "c-axis nodes" in underdoped 
and overdoped Ba(Fe-| - x Co x )2 A s 2 as suggested by c-axis 
heat transport experiments [381 ] . For the extremely hole- 
doped KFe 2 As2, the nodal superconductivity may result 
from the intraband pairing via antiferromagnetic fluc- 
tuations, due to the lack of electron pockets [l0| . For 
underdoped Bai_ a: K ;r Fe2As2, it is still not clear how 
the superconducting ga p transforms from nodeless to 



nodal at x ~ 0.16 [14(. The rest three compounds, 
BaFe2(Asi_a;Pa:)2, LaFePO, and LiFeP, have stimulated 
various interpretations of the effect of isovalent P substi- 
tution on the superconducting gap structure 22 - 25| . 

Our new finding of nodal superconductivity in 
Ba(Feo.64Ruo.36)2As2 reveals the similarity between the 
isovalently Ru- and P-substituted iron arsenides, there- 
fore the P substitution is not that special for inducing 
nodal superconductivity. In this sense, the mystery of 
P substitution in iron-based superconductors has been 
partially unwrapped. What next one needs to do is to 
find out whether there is a common origin for the nodal 
superconductivity in these isovalently substituted com- 
pounds. 

Due to the smaller size of P ion than As ion, one com- 
mon structural feature of the P-substituted compounds 
is the decrease of pnictogen height and increase of As-Fe- 
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As angle [32|, |4(j, |4l| . The substitution of larger Ru ion 
for Fe ion in Ba(Fei_a;Ru x )2As2 results in the increase 
of a lattice parameter and decrease of c lattice parame- 
ter, thus the decrease of pnictogen height and increase of 
As-Fe-As angle too [28}. Therefore, both the P and Ru 
substitutions cause the same trend of structure change in 
iron arsenides. 

LaFePO and LiFeP belong to the "1111" and "111" 
systems, respectively, and the P ions have fully substi- 
tuted As ions in LaFeAsO and LiFeAs. This is different 
from the partial P and Ru substitution in superconduct- 
ing BaFe2(Asi_ x P2;)2 and Ba(Fei_ x Ru x )2As2. In fact, 
both the fully substituted BaFe2P2 and BaRu2As2 are 
nonsuperconducting [ID, Another difference is the 
much lower T c of LaFePO and LiFeP, 6 and 4.5 K, re- 
spectively. For LaFePO, Kuroki et al. have attributed 
the low-T c nodal pairing to the lack of Fermi surface 7 
around (it, tt) in the unfolded Brillouin zone, due to the 
low pnictogen height [Hj]. Hashimoto et al. also related 
the nodal superconductivity in LiFeP to the pnictogen 
height [2l|. 

For Ba(Fei_ x Ru x )2As2 and BaFe2(Asi_ a; P; I ;)2, the 
substitutions start from the same parent compound 
BaFe2As2, and result in similar phase diagrams 2_| 29j | . 
The highest T c at optimal substitution, 20 and 30 K, are 
also close. Furthermore, the Fermi surface structures are 
roughly similar, with hole pockets around BZ center and 
electron pockets around BZ corners 17, [3(|. All these 



similarities suggest that the origin of the nodal supercon- 
ductivity in Ba(Fei_ x Ru :c )2As2 and BaFe_(Asi_ x P x )2 
may be the same. Suzuki et al. have proposed three- 
dimensional nodal structure in the largely warped hole 
Fermi surface and no nodes on the electron Fermi surface 
i25j. However, these seem inconsistent with the ARPES 
results, which have constrained the nodes on the electron 
pockets 17j. More careful considerations of the struc- 
tural parameters, band structure, and local interactions 
are needed to clarify whether there is a common origin 
for the nodal superconductivity in all these isovalently 
substituted iron arsenides. 

In summary, we have measured the thermal conductiv- 
ity of Ba(Fe .64Ruo.36)2As 2 single crystal down to 50 mK. 
A large kq/T at zero field and an H 1 / 2 field dependence 
of kq(H)/T at low field give strong evidences for nodal 
superconductivity in Ba(Feo.64Ruo.36)2As2. Comparing 
with previous P-substituted iron arsenides, our new find- 
ing suggest that the nodal superconductivity induced by 
isovalent substitutions may have the same origin, at least 
in Ba(Fei_ a; Ru a ;)2As2 and BaFe2(Asi_ K P a; )2. Finding 
out this origin will be important for getting a complete 
electronic pairing mechanism of the iron-based high-T c 
superconductors . 
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